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Abstract: 

Finite temperature lattice QCD indicates that the charmonium ground state J/'ip can 
survive in a quark-gluon plasma up to 1.5 or more, while the excited states Xc and 
xjj' are dissociated just above Tc. We assume that the Xc suffers the same form of sup¬ 
pression as that observed for the xjj' in SPS experiments, and that the directly produced 
J/ip is unaffected at presently available energy densities. This provides a parameter-free 
description oi J and suppression which agrees quite well with that observed in SPS 
and RHIC data. 


Recent studies of the behavior of charmonium states in a deconfined medium show that 
the ground state J/'0(1S) survives up to considerably higher temperatures than initially 
expected. In quenched QCD n-i, charmonium correlators show no signs of medium- 
induced suppression at least up to 1.5 Tc, while above 2 - 2.5 Tc, the signal is strongly 
modified or disappears. First work in QCD with two quark flavors supports these results 
[^. In contrast, the higher excited states seem to disappear very near T^; in quenched 
calculations, no signal for the Xc is seen at T = 1.1 Tc [S]. 

The results of direct spectral function studies are further supported by potential model 
analyses M-M, using the color-singlet free energy obtained in (quenched as well as un¬ 
quenched) lattice QCD to determine the heavy quark potential. These also lead to a 
J/'ijj dissociation temperature of 2 or higher, while Xc and xjj' disappear in the vicinity 
of 1.1 Tc. In contrast, earlier potential model work ca-iia. based on a heavy quark 
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interaction which nnderestimated the actnal QQ potential, had predicted a considerably 
lower J/il) dissociation temperatnre. 

Since J/'ip suppression was proposed as a signature for quark-gluon plasma formation in 
nuclear collisions PI, this modihcation of our understanding of the in-medium behavior 
of charmonia can be quite important for the interpretation of relativistic heavy ion data. 
Lattice calculations show that a temperature of 1.5 Tc corresponds to an energy density 
around 10 GeV/fm^, and 2 Tc to around 30 GeV/fm^, which could move the suppression 
of direct J/V’ production out of the range of RHIG. 

In hadron-hadron collisions m it is found that about 60% of the observed J/'ip's are 
directly produced as (IS) states, with the remainder coming to about 30% from Xc and 
10% from xp' decay. The hierarchy of suppression temperatures thus leads to a sequential 
suppression pattern cniis], with an early suppression of the xp' and Xc decay products 
and a much later one for the direct J/p) production. 

In this note, we want to consider the experimental results available now from the SPS and 
from RHIG, and show that the new theoretical understanding can be used to formulate a 
rather natural parameter-free description of the essential features of the data. 

Our considerations are based on the following scenario. The J/ip survival probability Sj/^ 
in A — A collisions is dehned as the ratio of the measured rate to that expected if the only 
modihcations are due to the presence of normal nuclear matter. We assume that Sj/^ 
consists of one term corresponding to the survival of directly produced J/tp^s and a 
second term Sx for those coming from the decay of the higher excited states Xc and pj', 

Sj/^ = 0.6 + OA Sx- (1) 

The relative contributions here are those observed in hadron-hadron collisions jHj. From 
the mentioned QGD studies we expect ~ 1 for energy densities up to 10 GeV/fm^ or 
more, while Sx is expected to show suppression effects around the deconhnement point, 
i.e., for e ~ 0.5 — 1.5 GeV/fm^. In principle, Sx could consist of two distinct terms, with 
different dissociation onsets for Xc and p>'. At present, however, neither calculational nor 
experimental accuracy seems to permit such hne structure studies, and we shall therefore 
combine the decay of the two states into one term. 

We hrst turn to the onset pattern of suppression and consider the SPS data for J/pj pro¬ 
duction from Pb—Pb and In—In interactions na. together with pj' data from Pb—Pb 
collisions UBI; the analysis of pj' production in the In—In data is not yet completed^. In 
addition, there are reference data from p — A collisions with several nuclear targets m, 
which dehne the necessary baseline for modihcations of the production due to normal 
nuclear matter. The combined effect of all possible modihcations was here parametrized 
in the form of nuclear absorption, leading to the absorption cross sections 

o'j/V' = 4.3 ± 0.3 mb (2) 

for the J/p) and 

= 7.1 ± 1.6 mb, (3) 

^We restrict ourselves here to symmetric {A—A) data and comment on the S — U results later on. 
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for the xjj', respectively ng. Using these in a Glauber analysis of A — A data provides 
the production rates {dai/dy)G, with i = J/t/’, ip', as they would be if there were no 
effects beyond those caused by the presence of normal nuclear matter m- The survival 
probability is then dehned as 

^ {dajdy) 

' {da,/dy)G^ ^ ’ 

describing whatever anomalous effects arise. The centrality dependence of the ^4 — ^4 
data is determined through the number Npart of participants, which is measured directly 
through a zero degree calorimeter. A Glauber analysis then provides the density Upart of 
participants in the transverse overlap region A of the collision pg, and the corresponding 
energy density is given by the Bjorken estimate 


_ / dNh \ _ UhWh 

Ato \ dy } To 


( 5 ) 


here {dNh/dy)AA denotes the hadron multiplicity at the given centrality, Wh the average 
hadron energy, and Vh the average number of hadrons emitted per participant nucleon (the 
values of and Wh can depend on centrality). For the equilibration time of the medium, 
we take Tq = 1 fm, so that corrections for other possible values can easily be carried out. 
In our context, however, the formation time of the charmonium states in question should 
be less than the formation time of the medium, which is the case if Tq = 1 fm. The actual 
values of e we will cite here were obtained by an event generator determination of the 
NA60 collaboration and is based on VENUS ini. It should be noted, however, that with 
constant ~ 2 and Wh — 0.5 GeV, we get very similar results, while an event generator 
determination based on RQMD as input (used by the NA50 collaboration |TH]) leads to 
values which are higher by about 10%. 

We now return to our basic scenario, assuming that at present energy densities the directly 
produced J/ip are unaffected, and the suppression patterns of the excited states Xc and 
Ip' are about the same. This implies that if we use the ip' data to form 0.4 + 0.6, then 

as function of the energy density this should coincide with the measured J/ip results. In 
Fig. m we see that this is indeed quite well fulhlled, for the overlap of J/ip and ip' data 
as well as for the convergence to the J/ip “saturation” value of about 60%. 



Figure 1: Universal ip' and J/ip suppression at the SPS 
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Next we want to check if this pattern continnes for higher energy densities and therefore 
tnrn to the recently presented preliminary RHIC data; its higher collision energy can 
provide correspondingly higher energy densities. The J/'ip prodnction rate Rau-Au in 
Au — Au interactions is given relative to the resnlt from scaled p—p collisions, as shown 
in Fig. 121 as fnnction of the nnmber of participant nncleons |21j . 
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Fignre 2: J/%Ij prodnction rates for Au —Au collisions at ^/s = 200 GeV j2I] 

In order to convert the rates Rau-Au into snrvival probabilities, we have to know what 
wonld be expected if only normal nnclear matter were present. At RHIC, this information 
is provided throngh d—Au stndies [ 22 ] ; the resniting nnclear modihcation factor, specifying 
the prodnction rate relative to scaled p—p collisions, is shown in Fig. |21 
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▼ Au-Au, RHIC, lyl=[1.2,2.2] Npart 


J_I_I_l_I_1_L 

100 200 300 


Fignre 3: J/'0 prodnction in d — Au collisions at y/s = 200 GeV 

To qnantify these RHIC resnits, with their presently rather limited statistics, we adopt 
a description similar to that nsed for SPS resnits and apply the well-known simplihed 
absorption form 

S ~ exp{-nocrdiss^}, (6) 

where L denotes the path of the cc in the nnclear medinm and hq = 0.17 fm“^ denotes 
normal nnclear density. A Glanber analysis [221 provides the relation between impact 
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parameter b and the number of collisions Ncou, and simple geometry gives L = — 

in terms of b and the nuclear radius Ra- A fit of Eq. (jU)) to the data of Fig. [Hlgives^ 

CTdissiy = 1-8) = 3.1 ± 0.2 mb 
<7diss{y = 0) = 1.2 ± 0.4 mb 

crdiss(2/ = -1.7) = -0.1±0.2mb (7) 


for the corresponding J/tp dissociation cross sections; for y = —1.7, there are thus es¬ 
sentially no nuclear modihcations. We note that here, as for the SPS case, these cross 
sections are just a global way to account for whatever nuclear effects can arise. A more 
detailed analysis based on shadowing and absorption is given in 1211 ; an analysis based on 
the Color Glass Condensate approach has recently been performed in I2nj- In the latter 
approach, the factorization of the shadowing and absorption corrections does not occur; 
nevertheless, here we use the equation just as a way to parameterize the data. 

For A —A collisions at RHIC energy, we make use of the same simplihed form (jHl). The 
geometry connecting the impact parameter b and path length L in p — Au and Au — Au 
collisions is illustrated in Fig. the relation between b and Npart is again given by a 
Glauber analysis ra. We thus here obtain for the survival probability 


sy(!/,iVp„,) 


7?j4a(?/5 A^part) 

exp{-no[crdiss(2/) + adiss(-|/)]T} ’ 


( 8 ) 


corresponding to the fact that for y ^ 0 the charmonium state passes one nucleus at 
rapidity y, the other at rapidity —y. 




Figure 4: Impact parameter relation between p —A and A —A collisions 

Applying eq. |Hlto the rates shown in Fig. [Altogether with the nuclear modihcation cross 
sections ca provides the survival probability as function of Npart- The corresponding 
energy densities have been calculated in a Glauber analysis based directly on the PHENIX 
Et data ra. and in Fig. El we compare the RHIG results to those from the SPS. 

It is seen that the two data sets are quite compatible, both in the onset and in the 
flattening at about 50 - 60%. Goncerning the RHIG data, it should be emphasized that 
the choice of tq = 1 fm is certainly debatable; a smaller value would move the RHIG 
points to correspondingly larger e values. 

^In the fit, we neglect the most peripheral point at N^oih which corresponds to 6 > Rau and is thus 
due to nuclear surface rather than medium effects. 
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Figure 5: J/'ijj suppression as function of energy density 


So far we have considered only symmetric {A —A) collisions. We find, however, that the 
xjj' production measured in S — U interactions at the SPS j2H] also agrees quite well with 
the pattern shown in Fig.^ In contrast, the reported S—U J/ip rates |2H1 do not show an 
onset of suppression at the centrality at which it sets in for In—In collisions. The reason 
for this is not clear, although two special features have been pointed out. The centrality 
dependence of the S—U data is determined by transverse energy {Et) measurements, not 
by the more reliable method based on the zero degree calorimeter specifying directly the 
number of spectator nucleons. For Pb — Pb collisions, it is observed that the centrality 
dependences obtained from and Ezdc measurements can in fact show differences. 
Moreover, it has been noted that a 10% shift in the normalization of the S—U data would 
lead to full agreement between all SPS data sets. 

J/^jj production at RHIC has also been adressed in terms of anomalous suppression fol¬ 
lowed by regeneration at hadronization [22] • Such a scenario assumes hrst a strong anoma¬ 
lous suppression of the overall J/tp production, including that of the IS* state, and subse¬ 
quently a renewed J/ip formation at the hadronization stage, due to a pairing of c andc 
quarks from different nucleon-nucleon collisions. The latter mechanism becomes possible 
at RHIC energies because of abundant cc production. It leads to rates increasing with 
centrality, which are taken to just compensate the dropping direct production. In such 
an approach, the agreement between central SPS data (with no regeneration) and RHIC 
rates (with considerable regeneration) is coincidental. We also note that the anomalous 
suppression assumed for direct J/pj production in the regeneration approach is not in 
accord with what we know today about J/ip survival in a quark-gluon plasma, as found 
in statistical QCD. 

Finally we turn to a further check of these considerations. It was pointed out some 
time ago that the effect oi Jjip suppression could also manifest itself in the transverse 
momentum behaviour EDI EH, and in fact the pattern resulting from sequential decay 
differs strongly from that due to regeneration inzi. 

The basic effect of a nuclear medium on the transverse momentum behaviour of hard 
processes is a collision broadening of the incident parton momentum; this in turn leads 
to a broadening of the transverse momentum distribution of the charmonia formed by 
hard parton interactions, (dominantly gluon fusion). It was shown that a random walk 
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approach leads to an average squared transverse J/V’ momentum 


{Pt)pA = {Pt)pp + 


( 9 ) 


ioT p —A and to 

{Pt)aa = {Pt)pp + (10) 

for y4 — y4 collisions. Here denotes the average number of pre-fusion collisions of 
the projeticle parton in the target nucleus A, and the sum of the average number 
of collisions of a projectile parton in the target and vice versa, at the given centrality. 
The parameter (5o specihes the average “kick” which the incident parton receives in each 
subsequent collision. The basic parameters determining the p^-broadening in nuclear 
matter are thus the elementary {p%)pp from p-p interactions and the value of determined 
by corresponding p — A data; both depend on the collision energy. The H-dependence of 
as well as the behaviour of as function of centrality can be obtained through a 
Glauber analysis; the latter dehnes the “normal” centrality dependence of {pt)aa- Such 
an analysis also has to include the normal absorption of the produced charmonia in 
nuclear matter; this effectively shifts the fusion point for the observed charmonia further 
“down-stream” EH- 

A compilation of J/'0 transverse momentum data from the SPS [SSI shown in Fig. 
ISl it clearly indicates hrst the increase of the average transverse momentum from p —p 
to p — A (for A = Pb), and then a further increase with centrality for nucleus-nucleus 
collisions. The preliminary data for the average p^ observed in J/'ip production at RHIC 
is shown in Fig. Here we note that while the muon arm data (|p| G [1.2, 2.2]) 

shows the expected broadening when going from p—p to d—Au, the central electron data 
(IpI < 0.35) does not follow this pattern. Since our analysis is based on such a broadening, 
we concentrate here on muon data. More statistics at central rapidity should clarify this 
problem. We note that since both {pt)pp and {pT)dA can in general depend on rapidity 
as well as on collision energy, each data set requires a separate analyis. 
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Figure 6: J/'if) transverse momentum 
behaviour at the SPS[H!?j 


Figure 7: J/il) transverse momentum 
behaviour at RHIC jsninii 


At SPS energy, one has {pt)pp = 1.25 ±0.05 (GeV/c)^ and {pt)pU = 1.49 ±0.05 (GeV/c)^ 
inn]- The average number of pre-fusion collisions is calculated in a Glauber analysis p?T| . 
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and with the normal nuclear absorption specihed by the average of eqns. (2/3) it is found 
to be about 3. From eq. o we then obtain 

5 SPS = 0.083 ± 0.023 GeV^ (11) 

for the average projectile parton broadening in the target nucleus. 

From the RHIC data, we obtain < >pp= 2.51 ± 0.21 (GeV/c)^ and < >dAu= 
3.96±0.28 (GeV/c)2 EHIMl; for the latter value, we have taken the average of the positive 
and negative rapidity ranges, since this is also done for the corresponding Au — Au data. 
A corresponding Glauber analysis, with normal nuclear absorption specihed by eq. ©, 
gives nearly 3.5 pre-fusion parton collisions and leads to 

^RHic ^ Q_42 ± 0.09 GeV^ (12) 

for the corresponding parton broadening in the large rapidity region. 

In a sequential dissociation scenario, the transverse momentum behaviour below the onset 
of exited state suppression is that of charmonia suffering only initial state broadening and 
normal nuclear absorption. Once the higher states are suppressed, one has once again 
only direct J/'ipA experiencing initial state effects and normal absorption. Hence apart 
from possible huctuations in the suppression region, one should observe the pt behaviour 
as given by eqs. © and m- In other words, the J/ip transverse momentum should be 
determined only by the initial nuclear medium. This again predicts a common behaviour 
of measurements from SPS and RHIG. Given the values of 60 as determined above, data 
for {p\) A A and {pt)pp dehne 

N^^ = {{pI)aa-{pI)pp}/5o (13) 

as a characteristic measure of transverse momentum behaviour. In Fig. |Hlwe show the 
SPS data from Pb—Pb |231 and In—In [23] collisions together with the RHIG muon data 
EH El and hnd that they indeed agree quite well. Once the corresponding broadening 
pattern for the RHIG electron data is determined, it should also follow this curve, even 
though the centrality dependent values for {p%)aa can be quite different. 

As mentioned, the centrality dependence of has also been calculated directly in a 
Glauber analysis EH; the result is included in Fig.|Hl It lies consistently somewhat higher 
than the results obtained from the data, which presumably comes from using a larger 
normal suppression. 

In contrast to the increasing p-r-broadening determined by initial state parton scattering, 
J/V’ production through cc pairing at hadronization leads to a centrality-independent 
{Pt)aa EH- Remnant direct production will of course modify this, but a strong regener¬ 
ation component should in any case considerably weaken the centrality dependence. 

The lack of the feed-down contributions to the observed J/'ijj and the presence of the 
plasma may also affect J/'ijj polarization [331137] . even though the theoretical description of 
quarkonium polarization has so far been notoriously difficult. Nevertheless, the predicted 
change of polarization may occur, and should be investigated experimentally. 
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Figure 8: Transverse momentum broadening at SPS and RHIC 

We conclude that present J/il) production data agree quite well with the expectations 
based on quark-gluon plasma formation. The observed onset of anomalous J/ijj suppres¬ 
sion now coincides, within errors, with that found for xjj' production, and the corresponding 
energy density agrees with that expected from finite temperature QCD for the dissoci¬ 
ation of higher excited charmonium states. The J/'ip production remaining beyond this 
initial anomalous suppression, at about 60 %, agrees with that predicted by a survival 
of directly produced IS* charmonium states and thus is also in accord with present QCD 
calculations. Further checks can come from measurements of tp' production in In —In 
collisions, from an eventual onset of direct J/p) suppression at higher e (LHC), and from 
similar results for T production in nuclear collisions. 
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